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ABSTRACT 


In  this  thesis,  the  adjacent  channel  interferer  ce  in  a  ultra  high  frequency  (UHF) 
satellite  channel  is  evaluated  by  simulation  and  differential  binary  phase-shift  keying 
(DBPSK)  is  compared  with  continous  phase  frequency-shift  keying  (CPFSK).  First,  a 
measure  of  the  interfering  power  is  obtained  and  a  method  to  compute 
carrier-to-interference  ratios  in  a  non-linear  channel  is  developed.  Next,  a  DBPSK 
receiver  is  simulated  when  two  interfering  channels  separated  in  frequency  are  present, 
and  bit  errors  are  detected  and  counted.  Then,  coherent  reception  of  minimum-shift 
keying  (MSK)  and  CPFSK  with  modulation  index  h=0.4  are  simulated  in  the  same 
conditions  as  DBPSK.  Finally,  noncoherent  MSK  is  analyzed  in  the  same  way  and  a 
comparative  behavior  is  obtained.  It  is  found  that  the  best  performance  in  the  presence  of 
adjacent  channel  interference  is  given  by  coherent  reception  of  MSK. 
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I.  INTRODUCTION 


A.  DISCUSSION 

The  main  goal  of  the  ultra  high  frequency  (UHF)  satellite  system  is  to  provide 
reliable  data  transmission  between  multiple  mobile  users.  In  a  digital  satellite  syston, 
performance  is  measured  in  terms  of  the  average  probability  of  bit  error.  Given  a 
sufficiently  large  bit-energy-to-single-sided-noise-power-spectral-density  ratio  (Eb/No), 
which  is  directly  proportional  to  the  carrier-to-noise  ratio  (C/N),  it  is  generally  assumed 
that  the  probability  of  bit  error  (Pb)  [Ref  1],  can  be  made  arbitrarily  small.  The  UHF 
satellite  is  a  frequency-division  multiple  access  (FDMA)  system.  Consequently,  when  a 
second  user  accesses  an  adjacent  chaimel,  some  spillover,  called  adjacent  channel 
interference,  will  occur,  and  this  will  degrade  the  performance  of  the  system,  even  for 
large  C/N,  since  the  effect  of  adjacent  channel  interference  is  to  reduce  C/N. 

a  INTERFERING  SOURCES  -  GENERAL  CONSIDERATIONS 

If  the  interference  source  is  assumed  to  be  a  statistically  independent  wide-sense 
stationary  random  process  of  zero  mean,  the  ovciall  carrier-to-noise-plus-interference 
ratio  can  be  expressed  by  [Ref  2], 

S=[(g)-‘+(7r'i-'.  (V 
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where  C/N  is  the  carrier-to-noise  ratio  of  the  overall  link,  and  C/I  is  the 
carrier-to-interference  ratio  of  the  overall  link  When  the  interferences  are  non-Gaussian 
but  numerous  and  none  of  them  has  a  dominant  effect,  their  joint  probability  density 
function  approaches  the  Gaussian  probability  density  function  as  stated  by  the  central 
limit  theorem.  The  effect  of  interference  in  this  case  can  therefore  be  assumed  to  be 
equivalent  to  the  effect  produced  by  a  single  additive  white  gaussian  noise  (AWGN) 
process  with  the  same  carrier-to-interference  ratio.  The  treatment  of  non-Gaussian 
interferences  as  equivalent  AWGN  generally  results  in  a  higher  predicted  probability  of 
bit  error  than  occurs  in  practice,  probably  because  the  sources  are  not  Guassian  and 
because  they  are  not  sufficiently  nume lous  for  the  central  limit  theorem  to  apply. 

The  consideration  of  interference  in  satellite  systems  is  of  utmost  importance.  The 
interference  could  come  from  such  different  sources  as  adjacent  satellite  systems, 
terrestrial  interference,  cross-polarization  interference,  adjacent  channel  interference,  and 
intermodulation  interference. 

Adjacent  satellite  system  interference  is  generated  by  an  earth  station  different  than 
the  one  rnder  consideration,  and  is  caused  by  the  power  received  through  the  antenna 
sidelobes  which  interferes  with  the  main  transmission.  This  effect  can  only  be  overcome 
by  designing  an  antenna  with  smaller  sidelobes. 

Terrestrial  interference  is  caused  by  terrestrial  networks  working  in  frequency 
bands  where  satellite  systems  have  channels  allocated.  In  the  case  of  the  UHF  satellite 
channel,  the  interference  could  come  from,  for  example,  terrestrial  mobile  systems  or 
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harbor  navigation  systems  It  is  known  that  these  kind  of  networks  have  a  limited  range, 
but  in  certain  conditions,  such  as  surface  ducts,  the  transmission  might  reach  unexpected 
distances  and  therefore  interfere  with  a  satellite  eartl>  station  that  is  located  outside  the 
area  of  influence  of  the  interfering  sources. 

Cross-polarization  interference  is  produced  in  satellite  systems  in  which  orthogonal 
linear  polarizations  are  employed  to  allow  frequency  reuse  The  depolarization  effect 
caused  by  rain  and  the  finite  cross-polarization  discrimination  of  the  earth  station  allow 
the  channels  to  interfere  with  one  another  in  spite  of  the  orthogonal  polarization  condition 
in  the  transmission  of  the  communication  message. 

Intermodulation  interference  is  caused  by  the  intermodulation  products  generated 
within  a  satellite  transponder  as  a  result  of  the  non-linear  amplification  of  multiple 
carriers  by  the  traveling  wave  tube  amplifier  (TWTA).  By  operating  the  high  power 
amplifiers  at  a  certain  output  backoff,  one  can  reduce  their  non-linear  effect  and  reduce 
the  intermodulation  interference. 

1.  Adjacent  Channel  Interference 

Another  source  of  interference  in  a  FDMA  satellite  link  is  the  adjacent  channel 
interference.  For  example,  the  power  spectral  density  of  binary  pnase-shift  keying 
(BPSK)  is  represented  in  Figure  1,  and  it  can  be  seen  that  most  of  the  energy  is 
concentrated  in  the  main  lobe  which  occupies  a  bandwidth  B  =  2/Tb  where  Tb  is  the  bit 
duration.  However,  the  sidelobes  of  the  spectrum  contain  some  energy  and  if  not  properly 
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filtered  out,  they  can  interfere  with  adjacent  channels  provided  the  separation  between 
them  is  not  high  enough.  This  situation  is  depicted  in  Figure  2. 

Obviously,  a  modulation  scheme  with  smaller  sidelobes  will  have  a  better 
performance,  as  far  as  adjacent  channel  interference  is  concerned,  than  one  with  higher 
sidelobes.  A  modulation  scheme  with  a  very  compact  mainlobe  and  low  sidelobes  is 
minimum-shift  keying  (MSK),  which  belongs  to  the  family  of  continuous  phase 
modulation  schemes  with  a  modulation  index  h=0.5.  The  basis  of  this  work  will  be  a 
comparative  analysis  of  the  adjacent  channel  interference  between  differential  binary 
phase-shift  keying  (DBPSK)  and  continuous  phase  frequency-shift  keying  (CPFSK),  a 
form  of  continuous  phase  modulation. 

2,  Jamming  Considerations 

The  interference  coming  from  a  jammer  can  be  considered  in  the  same  way  as 


interference  from  unintentional  sources.  That  is,  since 
c _ c_ 

K  “  N+J 

where  J=jamming  energy.  Then 

§  =  [(?)-■ +(7)-']-', 
where  C/J  is  the  carrier-to-jamming  ratio. 

Including  the  interference,  we  get 

g=[(Jr‘+(7r‘+(7)‘‘r‘- 


(2) 

(3) 


(4) 

The  term  (C/J)  '  is  called  the  jamming  margin  and  is  the  amount  of  jamming  the  system 
can  tolerate  for  a  certain  probability  of  bit  error.  Since  for  a  given  Pb  and  modulation 
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type,  a  unique  value  of  §  is  required,  it  also  determines  the  C/J  and  C/I  the  system  can 
accept  without  significantly  degrading  its  performance. 


Channel  1  Channel  2  Channel  3 

B  A  C 


f 

carrier 


Adjacent 

Channel 

Interference 

Figure  2.  Adjacent  Channel  Interference. 

C.  OBJECTIVE 

At  present,  the  UHF  satellite  described  in  the  Hughes  Aircraft  Company  Space  and 
Communication  Group  proposal  [Ref.  3]  cannot  successfully  be  used  at  bit  rates  of  4800 
bps  or  higher.  The  objective  of  this  thesis  is  to  demonstrate  that  this  limitation  is  due  to 
adjacent  channel  interference  and  can  be  solved  by  using  a  modulation  scheme  other  than 
DBPSK  such  as  MSK. 
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n.  ANALYSIS  OF  INTERFERING  POWER 


A.  UHF  SATELLITE  MODEL 

The  basic  model  of  the  satellite  channel  that  was  used  to  run  all  the  simulations 
contained  in  this  thesis  is  shown  in  Figure  3.  The  key  modules  were  adopted  from  [Ref 
3].  They  consist  of  a  prelimiter  filter,  a  hard  limiter  filter,  and  a  postlimiter  filter. 


Figure  3.  Basic  Sataiiita  Modai  with  Adjacent  Channei  Interference 


1.  Prelimiter  Filter 

The  prelimiter  filter  was  implemented  as  Chebyshev  Filter  with  6  poles  and  0.0 1 
dB  passband  ripple.  All  the  simulations  were  implemented  using  MATLAB.  MATLAB's 
filter  function  accepts  a  normalized  cutoff  frequency  value  between  0  and  1;  1 
corresponds  to  half  the  sampling  rate.  For  the  baseband  model  of  the  satellite  channel  a 
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sampling  frequency  (f.)  equal  to  384  kHz  was  chosen.  This  high  sampling  frequency  is 
required  since  f,  >  2  f,^  is  required  to  avoid  aliasing.  In  this  case  f,a„  =  100  kHz  is  the 
upper  frequency  of  the  upper  adjacent  channel  in  the  baseband  simulation. 

For  an  analog  cutoff  frequency  of  f^ ,  the  digital  cutoff  frequency  is 


f cutoff-  ■ 


(5) 


The  frequency  response,  both  magnitude  and  phase,  and  the  unit  impulse  response 
corresponding  to  this  filter  with  an  analog  cutoff  frequency  f^  =  12.57  kHz  are  plotted  in 
Figures  4,  5,  and  6,  respectively. 


Nonnatzad  Frequency 

Figura  5.  Phas*  Plot  of  Pralimttor. 


Tinw(eec)  -10* 

Figura  6.  Unit  Impulsa  Responaa  of  tha  Pralimiter  Fiiiar. 
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2.  Hard  Limiter 


The  hardlimiter  is  used  to  provide  constant  output  power  for  input  signal  power 
varying  from  the  noise  threshold  to  maximum  signal  input.  This  was  simulated  by 
dividing  each  sample  by  its  magnitude  such  that  each  complex  sample  is  on  the  unit 
circle. 

3.  Postlimiter  Filter 

This  filter  was  implemented  as  a  Chebyshev  filter,  with  4  poles  and  0.025  dB 
passband  ripple.  Based  on  the  same  considerations  as  before,  the  digital  cutoff  frequency 
for  this  filter  is  ~  0.0394,  since  the  analog  cutoff  frequency  f^  is  7.56  kHz.  The 
frequency  response,  both  magnitude  and  phase,  and  the  unit  impulse  response  for  this 
filter  are  plotted  in  Figures  7,  8,  and  9,  respectively. 


Nomnfaad  Fraquancy 

Figuro  7.  Froquency  Rosponso  of  the  Postlimitor  Filter. 
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Normiizad  Fraquancy 

Figur*  8.  Phass  Plot  of  tho  Postlimitor  Filter. 


'nm*(Mc)  10* 

Figure  9.  Unit  Impulse  Response  of  the  Postlimiter  Filter. 
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B.  LINEAR  ESTIMATION  OF  INTERFERING  POWER 


Having  described  the  basic  components  of  the  UHF  satellite  channel  we  are  now  in 
a  position  to  analyze  the  interference  from  adjacent  channels.  Initially,  only  interference 
due  to  the  upper  channel  is  considered.  Consequently,  consider  the  designated  channel  to 
be  a  baseband  channel  with  C«n«r~  0  (see  Figure  2).  The  results  can  be  easily  extended  to 
more  than  one  channel. 

The  separation  in  frequency  between  channels  plays  an  important  role.  Not  all  the 
channels  of  the  UHF  satellite  are  equally  spaced  in  frequency.  The  worst  case,  a 
frequency  separation  equal  to  100  kHz,  was  used  in  the  simulation.  The  first  experiments 
used  DBPSK  as  the  modulation  scheme.  A  block  diagram  for  the  experiment  is  shown  in 
Figure  10. 


Baseband 

Channel 


Branch  1 


Branch  2 


Adjacent 

Channel 


Figure  10.  Block  Diagram  of  Simulation. 


From  Branch  1,  the  power  in  the  baseband  from  the  on-charmel  signal  was 
computed.  Similarly  from  Branch  2,  the  power  in  the  baseband  coming  from  the  adjacent 
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channel  was  obtained.  The  results  for  this  experiment  are  presented  in  Table  I  This 
experiment  gives  a  first  indication  of  the  interfering  effect,  but  it  is  not  useful  to  provide 
an  accurate  value  of  the  C/I  ratio.  No  consideration  was  given  to  the  correlation  of  the 
processes  introduced  by  the  limiter  since  both  channels  were  analyzed  separately. 


TABLE  1.  NORMALIZED  BASEBAND  AND  ADJACENT  POWER 

FOR  DBPSK 


2400  bps 

4800  bps 

Bit  Rate 

9600  bps 

19200  bps 

ADJACENT 

POWER 

0.11 

0.19 

0.33 

0.43 

BASEBAND 

POWER 

0.97 

0.94 

0.89 

0.77 

BP/AP 

(dB) 

9.37 

6.96 

4.31 

2.56 

A  brief  look  at  Table  1  shows  that  the  figures  obtained  are  as  expected.  For  a  higher 
bit  rate  the  power  spectral  density  of  DBPSK  is  wider  [Ref  2];  more  power  from  the 
adjacent  channel  and  less  of  the  baseband  power  is  in  the  baseband  channel  bandwidth. 

For  the  second  experiment,  continuous  phase  modulation  was  selected  as  a  possible 
scheme  for  improvement  with  regard  to  adjacent  channel  interference.  The  same 
simulation  was  run,  and  the  results  for  MSK  (CPFSK  with  a  modulation  index  h  =1/2) 
are  shown  in  Table  2. 
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TABLE  2.  NORMALIZED  BASEBAND  AND  ADJACENT  POWER 

FOR  MSK 


2400  bps 

Bit  Rate 

4800  bps  9600  bps 

19200  bps 

ADJACENT 

POWER 

0.09 

0.15 

0.27 

0.36 

BASEBAND 

POWER 

0.99 

0.99 

0.99 

0.97 

BP/AP 

(dB) 

10.67 

8.19 

5.67 

4.27 

Two  other  attempts  were  made  to  find  out  if  a  different  modulation  index  h  could 
improve  performance.  CPFSK,  with  indexes  ranging  from  0.1  to  1,  was  analyzed  and  the 
results  are  plotted  in  Figures  11  and  12.  Similar  performance  is  expected  for  MSK  and 
CPFSK  with  h  =  0.4.  However,  a  small  improvement  can  be  detected  at  19200  bps  for 
h=0.4.  Therefore,  the  simulation  was  run  for  CPFSK  with  h  =  0.4  and  the  results  can  be 
seen  in  Table  3. 
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TABLE  3.  NORMALIZED  BASEBAND  AND  ADJACENT  POWER 

FOR  CPFSKWITH  h-0.4 


2400  bps 

Bit  Rate 

4800  bps  9600  bps 

19200  bps 

ADJACENT 

POWER 

0.09 

0.15 

0.27 

0.36 

BASEBAND 

POWER 

0.99 

0.99 

0.99 

0.98 

BP/AP 

(dB) 

10.56 

8.19 

5.58 

4.35 

ModuWkin  Indn  h 

Figurall.  BasAband  Powsr. 
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Moduiltion  mdax  h 


Figur«12.  Adjacent  Powar. _ _ 

The  final  experiment  used  Gaussian  MSK  as  a  modulation  scheme.  This  particular 
type  of  modulation  is  fully  described  by  Murota  and  Hirade  [Ref  4],  It  is  stated  to  have  a 
better  performance  than  MSK  in  certain  aspects,  such  as  ISI  degradation.  The  simulauon 
was  therefore  run  for  this  particular  case,  and  the  results  are  shown  in  Table  4. 
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TABLE  4.  NORMALIZED  BASEBAND  AND  ADJACENT  POWER 

FOR  GAUSSIAN  MSK 


2400  bps 

Bit  Rate 

4800  bps  9600  bps 

19200  bps 

ADJACENT 

POWER 

0.07 

0.17 

0.31 

0.55 

BASEBAND 

POWER 

0.99 

0.99 

0.99 

0.97 

BP/AP 

(dB) 

11.33 

7.73 

5.02 

2.47 

Comparing  all  the  results  obtained  so  far,  it  can  be  concluded  that  MSK  and  CPFSK 
with  h  =  0.4  are  candidates  to  outperform  DBPSK  in  the  rase  of  adjacent  channel 
interference.  Therefore,  a  more  detailed  study  is  necessary  to  obtain  a  more  accurate 
estimate  of  the  actual  carrier-to-interference  ratio.  An  approach  to  deal  with  this  situation 
is  developed  in  the  next  section. 


C.  CARRIER-TO-INTERFERENCE  RATIO  FOR  A  NON-LINEAR  CHANNEL 
Because  of  the  presence  of  the  hard  limiter  in  the  satellite,  the  system  is  not  linear, 
and  therefore  a  more  accurate  technique  to  estimate  the  carrier-to-interference  ratio  is 
necessary.  The  method  chosen  consists  of  estimating  the  on-channel  signal  and  removing 
it  from  the  on-channel  plus  interference  signals  in  order  to  estimate  the  interference.  The 
block  diagram  in  Figure  13  illustrates  this  technique. 
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Figure  13.  Block  Diagram  of  Estimation  Method. 


Since  the  adjacent  channel  signal  is  being  generated  independently  of  the 
on-channel  signal,  the  input  processes  are  uncorrelated  with  one  another.  The  best 
estimate,  Y,,  of  the  on-channel  signal  in  Y^occurs  when  Y,  is  orthogonal  to  the  error, 
Ya-Yj.  This  is  when 

EWA-Yi)xYi\  =  0  ,  (6) 

which  leads  to 


miYA-X,xG)xXiG\  =  0  . 


(7) 


and 


G  = 


E(.Ya*X,) 

EPCJ) 


(8) 


The  carrier-to-interference  power  ratio  can  then  be  expressed  as 
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This  technique  was  used  to  estimate  the  carrier-to-interference  ratio  for  DBPSK,  MSK, 
and  CPFSK  with  h  =  0.4.  The  results  are  presented  in  Table  5. 


TABLE  5.  CARRIER  TO  INTERFERENCE  RATIOS  (IN  dB) 
_ FOR  DBPSK,  MSK,  AND  CPFSK _ 

Bit  Rate 


Mod.  Scheme 

2400  bps 

4800  bps 

9600  bps 

19200  bps 

DBPSK 

25.5 

22.66 

18.14 

16.24 

MSK 

30.48 

27.26 

24.48 

19.89 

CPFSK 

(h=0.4) 

30.58 

27.43 

24.86 

21.03 

In  the  same  way,  and  based  on  the  independence  assumption  among  channels,  the 
carrier-to-interference  ratio  for  two  adjacent  channels  can  be  calculated.  Table  6  shows 
the  results  for  DBPSK  and  CPFSK  (h  =  0.4). 
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TABLE  6.  CARRIER  TO  INTERFERENCE  RATIOS  (IN  dB) 

_ FOR  DBPSK  AND  CPFSK  -  TWO  ADJACENT  CHANNELS _ 

Bit  Kate 

Mod.  Scheme  2400  bps  4,800  bps  9600  bps  19200  bps 
DBPSK  22.23  19.09  14.89  12.46 

CPFSK  27.58  24.43  21.78  17.98 

(h  =  0.4) 

The  results  in  Table  6  were  obtained  by  locating  a  lower  interference  channel  100 
kHz  from  the  baseband  channel.  The  simulation  was  then  run  with  both  upper  and  lower 
interfering  channels. 

These  results  show  that  either  MSK  or  CPFSK  (h  =  0.4)  have  C/I  significantly 
higher  than  DBPSK  at  all  data  rates.  The  procedure  can  be  continued  by  adding 
additional  channels  spaced  in  frequency  by  100  kHz  from  the  on'Channel  signal. 
However,  it  is  assumed  that  the  total  adjacent  channel  interfering  power  is  dominated  by 


the  first  adjacent  channels. 


m.  ANALYSIS  OF  THE  SATELLITE  CHANNEL  FOR  DIFFERENT 
MODULATION  TECHNIQUES 

A.  DBPSK  ANALYSIS  AND  SIMULATION  RESULTS 

For  this  modulation  scheme,  a  model  similar  to  the  one  used  by  Khanaman  [Ref.  S] 
Avas  simulated.  First,  a  lower  and  upper  interfering  channel  separated  100  kHz  in 
frequency  from  the  baseband  channel  were  simulated.  Since  the  computed 
carrier-to-interference  ratio  (C/I)  for  this  case  is  very  high  (see  Table  5),  no  errors  were 
expected  to  be  found  due  to  the  adjacetit  channels.  The  limitations  imposed  by  the 
computer  simulation  run  time  (no  more  than  1000  bits  were  simulated)  do  not  allow  the 
channel  to  be  analyzed  in  the  region  where  the  probability  of  bit  error  is  expected  to  be  as 
low  as  10'*^.  Therefore,  it  was  decided  to  reduce  the  frequency  separation  so  as  to  cause 
some  errors  to  appear  in  order  to  have  a  measure  to  compare  DBPSK  and  CPFSK. 

Obviously,  the  count  of  the  number  of  errors  in  any  Monte  Carlo  simulation  does 
not  represent  accurately  the  probability  of  bit  error,  because  only  a  finite  number  of  trials 
are  possible.  However,  the  number  of  errors  can  provide  a  good  idea  of  comparative 
behaviour  between  two  different  modulations  when  the  same  parameters  are  used  for  the 
channels. 

Consequently,  a  second  simultation  was  run  placing  two  adjacent  channels  at  +/-  IS 
kHz  and  a  third  simulation  was  run  locating  the  interfering  sources  at  +/-  12.5  kHz.  To 
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have  even  more  data  to  analyze,  the  channel  was  tested  for  three  different  Eb/No 
conditions;  14  dB,  12  dB,  and  10  dB  (in  the  last  one  +/>  25  kHz  was  used  instead  of  +/- 
15  kHz). 

A  block  diagram  of  the  channel  is  presented  in  Figure  14,  and  the  simulation  results 
are  presented  in  Tables  7,  8,  and  9.  The  different  codes  that  were  used  to  simulate 
DBPSK  can  be  found  in  Appendix  A. 


Figure  14.  Block  Diagram  of  DBPSK  Satallita  Channel  Simulation. 


TABLE  7.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  WITH 


14  dB  Eb/No. 

2400  bps 

4800  bps 

9600  bps 

19200  bps 

ADJACENT 
CHANNELS 
(100  kHz) 

0 

0 

0 

0 

ADJACENT 
CHANNELS 
(15  kHz) 

0 

0 

0 

90 

ADJACENT 
CHANNELS 
(12.5  kHz) 

0 

0 

1 

147 

TABLE  8.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  WITH 


12  dB  Eb/No. 

2400  bps 

4800  bps 

9600  bps 

19200  bps 

ADJACENT 
CHANNELS 
(100  kHz) 

0 

0 

0 

2 

ADJACENT 
CHANNELS 
(15  kHz) 

0 

0 

2 

100 

ADJACENT 
CHANNELS 
(12.5  kHz) 

0 

1 

4 

148 
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TABLE  9.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  WITH 


10  dB  Eb/No. 

2400  bps 

4800  bps 

9600  bps 

19200  bps 

ADJACENT 
CHANNELS 
(100  kHz) 

0 

0 

0 

3 

ADJACENT 
CHANNELS 
(25  kHz) 

0 

0 

0 

12 

ADJACENT 
CHANNELS 
(12.5  kHz) 

1 

1 

8 

152 

B.  MSK  ANALYSIS  AND  SIMULATION  RESULTS 
1.  Coherent  Reception 

In  coherent  MSK,  it  is  assumed  that  the  initial  phase  of  the  transmitted  signal  is 
perfectly  known  at  the  receiver.  Two  basic  coherent  receivers  were  modeled  for  this 
study.  The  first  is  explained  by  Haykin  [Ref  6].  Essentially,  it  consists  of  a  correlator 
receiver  with  two  branches  where  the  decision  is  made  by  alternatively  evaluating  the 
signal  after  integrating  it  over  a  period  equal  to  twice  the  bit  duration  (2  Tb)  with  one  bit 
offset.  A  simplified  block  diagram  can  be  seen  in  Figure  15.  The  decision  logic  is  shown 
in  Figure  16. 


Figura  15.  Coharant  Damodulator  of  MSK. 

I 

I 


X2 


Xi 


♦ 


Figura  16.  Dacision  Logic  for  MSK  Racaivar. 
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The  same  procedure  was  followed  with  DBPSK.  The  simulation  was  run  placing 
the  adjacent  channels  at  +/-  100,  +/-  15,  and  +/-  12.5  kHz.  (As  before,  when  a  10  dB 
Eb/No  was  used,  the  frequency  spacing  was  +/-  25  kHz  instead  of  +/-  15  kHz.)  The 
results  can  be  seen  in  Tables  10,  1 1,  and  12.  It  follows  that  there  is  an  improvement  in  the 
system  if  MSK  is  used  because  no  errors  were  found  until  the  channels  were  unacceptably 
close,  and  even  in  this  situation  the  number  of  errors  computed  was  considerably  lower 
than  in  the  case  of  DBPSK. 


TABLE  10.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  WITH 

14  dB  Eb/No. 


ADJACENT 
CHANNELS 
(100  kHz) 

ADJACENT 
CHANNELS 
(15  kHz) 

ADJACENT 
CHANNELS 
(12.5  kHz) 


2400  bps  4800  bps  96CX}bps  19200  bps 
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TABLE  11.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  WITH 

12  dB  Eb/No. 


ADJACENT 
CHANNELS 
(100  kHz) 

ADJACENT 
CHANNELS 
(15  kHz) 

ADJACENT 
CHANNELS 
(12.5  kHz) 


2400  bps  4800  bps 


9600  bps 
0 

0 

0 


19200  bps 
0 


6 


71 


TABLE  12.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  WITH 

10  dB  Eb/No. 


ADJACENT 
CHANNELS 
(100  kHz) 

ADJACENT 
CHANNELS 
(25  kHz) 

ADJACENT 
CHANNELS 
(12.5  kHz) 


In  the  second  receiver,  the  Viterbi  algorithm  is  used  to  decode  the  MSK  signal. 
As  explained  by  Proakis  [Ref.  7],  the  states  for  the  phase  of  MSK  can  be  +/-  7t/2,  0,  and  tc. 
The  number  of  states  can  be  reduce  if  the  signal  is  premultiplied  by  e^".  Note  that  an 
MSK  signal  leaving  from  a  phase  of  zero  will  increase  the  phase  hyidl  if  the  input  is  a 
logical  "one."  The  effect  of  the  premultiplier  adds  another 7t/2,  which  leads  to  a  fmal 
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phase  of  n.  If  the  input  is  a  logical  “zero",  the  phase  will  decrease  by  -idl.  The  effect  of 
the  premultiplier  leads  to  a  fmal  phase  state  of  zero.  In  other  words,  the  premultiplier 
reduces  the  number  of  phase  states  from  four  (+/-  idl,  0, 7t)  to  two  (0,  n).  The  trellis 
phase  diagram  for  MSK  with  premultiplication  is  illustrated  in  Figure  17. 


State 


2  input=*0"  2  ’O" 


Figure  17.  Phase  Trellis  Diagram  for  MSK  After  Premultipiication. 

After  the  premultiplication,  the  received  signal  is  correlated  and  each  branch  is 
used  as  an  input  to  the  Viterbi  algorithm,  as  illustrated  in  the  receiver  block  diagram 
shown  in  Figure  18.  A  soft  Viterbi  algorithm  tracks  the  phase  changes  along  the  trellis 
and  decides  on  the  most  probable  path  by  considering  as  a  decision  rule  the  minimum 
euclidean  distance  to  the  four  points  in  the  two-dimensional  (2-D)  plane  formed  from  the 
receiver  output  pairs.  In  this  receiver  four  output  pairs  are  possible,  depending  on  the 
previous  phase  state  and  the  input  bit  (see  Figure  1 9). 
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exp(/(27t/‘20) 


Figura  18.  R«c«iv«r  Block  Diagram  for  VKarbi  Docoding. 


Without  noise,  the  output  pair  will  coincide  exactly  with  one  of  the  four  possible 


points,  depending  on  the  input  bit  and  the  previous  phase  state,  as  described  in  Table  13. 
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l  ABLE  13.  POSSIBLE  OUTPUT  POINTS  FROM  THE 
VITERBI  DEMODULATOR 


Previous  State 


It  is  useful  to  notice  that  if  the  input  is  a  "I”,  independent  of  the  previous  phase 
state,  the  output  pair  will  be  on  the  "y"  axis,  whereas  if  the  input  is  a  "0"  the  output  will  be 
on  the  "x"  axis. 

The  simulation  was  run  using  MSK  and  the  Viterbi  receiver.  The  results  were 
found  to  be  a  bit  degraded  (-0.5  dB)  with  respect  to  the  receiver  described  in  Figure  16, 
but  still  superior  to  DBPSK.  The  different  codes  that  were  used  to  simulate  MSK  can  be 
found  in  Appendix  B. 

2.  Coherent  Reception  of  CPFSK  with  h=0.4 

It  was  seen  in  the  previous  chapter  that  CPFSK  with  modulation  index  h=0.4 
increases  the  C/I  by  a  small  amount  and  could  therefore  lead  to  better  performance  as  far 
as  this  interference  is  concerned. 

Since  CPFSK  with  h=0.4  is  not  an  orthogonal  signaling  set  [Ref  1],  the  first 
coherent  receiver  that  was  used  to  decode  MSK  cannot  be  used  as  a  demodulator. 
However,  the  Viterbi  algorithm  can  still  determine  a  maximum  likelihood  path  through 
the  phase  trellis  diagram  and  optimally  decode  the  signal.  A  Viterbi  receiver  was 
designed  for  the  h=0.4  CPFSK  signal.  For  this  signal,  there  are  five  possible  phase  states 
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(0,  +/-  27t/5,  +/-  4n/5),  and  th  ..mot  be  reduced  by  premultiplication  The  possible 
output  pairs  in  the  2-D  euclidean  plane  are  therefore  10.  The  large  number  of  points  leads 
to  a  serious  degradation  of  receiver  performance  since  the  points  on  the  euclidean  plane 
are  very  close  to  one  another.  When  noise  is  added,  a  very  high  signal-to-noise  ratio  is 
required  to  avoid  performance  degradation.  Since  this  is  not  the  case  for  a  satellite 
channel,  CPFSK  with  h=0.4  cannot  perform  as  well  as  MSK  even  though  it  has  a  very 
small  advantage  with  respect  to  adjacent  channel  interference.  The  code  that  was  written 
to  simulate  CPFSK  (h=0.4)  can  be  found  in  Appendix  C. 

3.  Noncoherent  Reception  of  MSK 

Coherent  reception  is  difficult  to  carry  cut  in  terms  of  receiver  complexity 
because  carrier  synchronization  is  required.  It  was  decided  to  investigate  the  performance 
when  noncoherent  MSK  is  used.  Several  noncoherent  receivers  have  been  described  in 
the  literature  [Ref  9,  10,  11].  The  best  performance  against  noise  is  obtained  by  using 
the  noncoherent  receiver  developed  by  Crozier,  et.  al.  [Ref  11].  A  block  diagram  of  the 
receiver  can  be  seen  in  Figure  20. 
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The  receiver  consists  of  a  differential  detection  branch  that  measures  the 


difference  in  phase  between  two  successive  signaling  intervals,  and  a  second  branch 


where  the  symbol  detected  from  the  difference  in  phase  between  two  alternate  signaling 


intervals  can  be  interpreted  as  the  parity  check  sum  of  two  successive  transmitted  data 
elements.  These  two  symbols  correspond  to  data  and  parity  of  a  rate  1/2 
single-error-correcting  self-orthogonal  convolutional  code;  therefore,  performance  can  be 
improved  by  using  the  decoder  for  this  error  correcting  code  [Ref.  9], 

To  get  even  better  performance,  two  filters  are  added.  The  reception  filter  is  a 
4-pole  phase  equalized  Butterworth  filter  with  filter-bandwidth-bit-duration  product 
(BT)=1.1  and  the  demodulation  filter  is  a  4-pole  phase  equalized  Butterworth  filter  with 
BT=1.5. 

The  results  of  the  simulation  can  be  seen  in  Tables  13,  14,  and  IS. 
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TABLE  13.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  VyATH 

14  dB  Eb/No. 


ADJACENT 
CHANNELS 
(100  kHz) 

ADJACENT 
CHANNELS 
(15  kHz) 

ADJACENT 
CHANNELS 
(12.5  kHz) 


2400  bps 

0 


0 


1 


4800  bps 

0 

0 

2 


9600  bps 

0 

0 

5 


19200  bps 

1 

84 

157 


TABLE  14.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  WITH 

12  dB  Eb/No. 


ADJACENT 
CHANNELS 
(100  kHz) 

ADJACENT 
CHANNELS 
(15  kHz) 

ADJACENT 
CHANNELS 
(12.5  kHz) 
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TABLE  15.  NUMBER  OF  ERRORS  FOR  THE  SATELLITE  CHANNEL  WITH 


10  dB  Eb/No. 

2400  bps 

4800  bps 

9600  bps 

1 9200  bps 

ADJACENT 
CHANNELS 
(100  kHz) 

0 

0 

0 

1 

ADJACENT 
CHANNELS 
(25  kHz) 

0 

0 

0 

2 

ADJACENT 
CHANNELS 
(12.5  kHz) 

3 

10 

20 

174 

From  these  results,  it  can  be  concluded  that  the  noncoherent  receiver  would  work 
in  high  signal-to-noise  ratio  situations.  However,  in  a  noisy  channel  the  single  error 
correction  circuit  cannot  correct  the  data  transmitted,  and  the  receiver  cannot  perform 
even  as  well  as  DBPSK. 

The  code  written  to  simulate  noncoherent  MSK  can  be  found  in  Appendix  D. 

4.  Coherent  MSK  Revisited 

The  coherent  reception  of  MSK  needs  both  a  carrier  recovery  circuit  and  a  clock 
recovery  circuit.  An  example  of  a  circuit  suitable  for  this  purpose  can  be  found  in  the 
work  of  deBuda  [Ref  12].  The  original  circuit  generates  90  degree  phase  and  multiples 
of  90  degree  phase  ambiguity  in  the  reference  canier  phase.  An  improvement  to  the 
circuit  that  resolves  the  phase  ambiguity  of  +/-  90  degrees  can  also  be  found  in  deBuda's 
work.  One  way  to  solve  the  remaining  180  degree  phase  ambiguity  is  by  differentially 
encoding  the  bit  stream.  However,  this  last  step  is  not  necessary  since  the  Viterbi 


algorithm  resolves  this  ambiguity  automatically.  The  trellis  diagram  remains  the  same 
when  this  ambiguity  is  introduced  in  the  receiver,  but  the  0  and  n  phase  states  are 
interchanged,  as  shown  in  Figure  21 . 


Assume  the  coherent  references  are  shifted  incorrectly  by  180  degrees.  The 


situation  is  pictured  in  Figure  22. 


exp(/(2ji/2/+7t)) 


Figure  22.  Viterbi  Demodulator  with  Coherent  Referencee  Shifted  by  180  Degrees. 


The  output  of  each  branch  will  be: 

jr=  Re[s^{t)  *  expo  *  (27t/ir + x))P  ;  [10J 

X=\^Re[m+jQiO)*(cos(2itfit+n)+jsm(2itfit+n))]a ;  [11] 

X=\^m<^s(2irfxt+n)-Q(t)sm(2itftt+n))dt ;  [12] 

-[Jo'"(/(0cos(2j5f,/)  -  Q(t)sm(2itfii))dt]  [13] 

In  the  same  way, 

Y=-\i^mcos(2nf2t)-Q0)sin(2iif2t))dt]  .  [14] 


But  the  terms  in  brackets  are  the  outputs  of  the  demodulator  if  the  phases  are  not  shifted. 
Therefore,  since  X  and  Y  are  the  components  in  the  2-D  plane  of  the  output  point,  it  is 
easy  to  see  that  the  new  output  has  been  shifted  by  180  degrees.  From  Table  13,  if  the 
ou^ut  is  shifted  by  180  degrees,  the  Viterbi  algorithm  still  decodes  it  as  the  same  bit. 
Only  180  degree  ambiguities  can  be  resolved  in  this  fashion. 

In  summary,  the  clock  and  carrier  recovery  circuit  can  be  implemented  as  shown 
by  deBuda  [Ref  12],  however  it  is  not  necessary  to  differentially  encode  the  message  if  a 
Viterbi  algorithm  is  used  as  the  decoder. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  results  obtained  in  this  thesis  show  that  coherent  minimum-shift  keying  with 
Viterbi  decoding  can  improve  the  performance  of  a  UHF  satellite  system  when 
interference  coming  from  adjacent  channels  is  the  main  concern. 

It  was  shown  that  continuous  phase  frequency-shift  keying  with  modulation  index 
other  than  h=0.5  and  non-coherent  reception  of  MSK  are  not  suitable  since  in  one  way  or 
another  their  performance  is  seriously  degraded  in  a  noisy  environment. 

For  coherent  MSK,  a  carrier  recovery  circuit  that  does  not  add  great  complexity  to 
the  receiver  and  that  does  not  adversely  affect  the  performance  of  the  coherent  MSK 
modulation  is  required.  A  circuit  was  presented  that  satisfies  these  criteria.  It  was 
demonstrated  that  the  circuit's  residual  180  degree  phase  ambiguity  is  solved  by  the 
Viterbi  algorithm  without  differentially  encoding  the  data. 

Unfortunately,  the  results  obtained  in  this  work  do  not  fully  support  the  thesis  that 
adjacent  channel  interference  is  limiting  satellite  channel  bit  rate  since  no 
interference-caused  errors  are  observed  in  the  simulation  when  the  chaimels  are  separated 
by  100  kHz.  The  results  are  consistent  v  ith  the  limitations  of  the  computer  model  and  the 
high  carrier-to-noise  ratios  computed  for  this  system,  and  the  work  done  for  this  thesis 
provides  a  good  comparative  idea  of  the  behavior  of  the  channel  under  those 
circumstances.  However,  it  is  necessary  to  have  a  more  accurate  tool  to  measure  the 
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actual  performance  of  the  satellite  to  determine  whether  coherent  MSK  has,  in  fact,  any 
real  benefits  for  UHF  satellite  communications. 
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APPENDIX  A. 


AWGN.M  CAWGN  FUNCTION) _ 

function  y  «  awgn(x, sigma) 

%  Awgn  is  an  M_file  that  adds  awgn  to  the  matrix  x.  The  standard  deviation  of 
%  the  noise  is  also  an  input  (sigma)  and  it  has  to  be  change  according  to  the 
%  different  Eb/No  that  are  desired  to  simulate,  v^«^e  Eb/No  »  l/(2*sigma^2). 
Irr.ccJ »  sizefx): 
seed  s  0; 
rand(’nonnai’); 
rand('seed’,seed); 
w  =  rand(rr,cc)  +  j*rand(rr,cc); 
y  «  X  +  sigma.  *w; 


COMPARE.M  (NUMBER  OF  ERRORS  FUNCTION) _ 

function  out  >  compare(in,in1) 

%  This  M_file  accepts  two  vectors  of  equal  length  composed  by  zeros  and  ones 
%  and  returns  the  number  of  bits  in  which  both  vectors  do  not  agree, 
com  s  abs(in  •  in1); 
out  s  sum(com); 


DBPSK.M  (DBPSK  MAIN  PROGRAM) _ 

%  receiver  for  DBPSK 
m«  1002; 

md_o1  s  msg(40,m);  %  Creating  the  random  message 
md_o2  s  msg(43,m);  %  Creating  the  interference  sources. 
rrKl_o3  s  msg(65,m); 

dif.ol  s  dif_cod(md_o1);  %  Differentially  encocHng  the  message. 
dif_o2  «  dif_cod(md_o2);  %  Differentially  encoding  the  interfering 

%  messages 

dlf_o3  «  dlf_ood(md_o3): 

map_o1  s  map(dif_o1);  %  Mapping  the  message 
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map_o2  «  map(dif_o2);  %  Mapping  the  inteffering  message 

map_o3  s  map(dif_o3): 

dd  a  [38  38  38  41];  %  Filter  delays 

T  a  [(1/2400)  (1/4800)  (1/9600)  (1/19200)];  %  Bit  durations 

t  a  1/384000;  %  Sampling  interval. 

f1  »0; 

deltaj  a  100000;  %  Frequency  separ^ion. 
bit  a  [160  80  40  20]; 

sigma  a  [2*s(^(2)  2  8qrt(2)  1];  %  Standard  deviation  of  the  noise 
clear  dif_o1  dlf_o2  dif_o3  md_o2  rrKf_o3 
for  ja4:4. 

mod_sig1  a  modul(map_o1,T(j).t,f1);  %  BPSK  modulation 

mod_sig2  a  modul(map_o2,T0‘).t.dolta_f): 

mod_sig3  a  modui(map_o3,T0').t.-delta_f); 

mod_sig  a  mod.sigi  +  mod_sig2  +  mod_sig3;  %  Adding  the  signals 

clear  mod_sig1  mod_sig2  mod_8ig3 

ch_sig  a  awgn(mod_sig,sigma(j));  %  Adding  the  Gaussian  noise 

clear  mod_sig 

ch_sig  a  ch_sig’; 

ch_sig  a  ch_8ig(:); 

ch_8lg  a  ch_8lg'; 

[b1,a1]  a  cheby1(6,.01,0.0651); 

prefiLsig  »  filter(b1  ,a1  ,ch_sig);  %  Prefiltering  the  signal 

clear  ch.sig 

iim_8ig  a  limiter(prefiLsig);  %  Hard  limiter  effect 
clear  prefiLsig 

[b2,a2]  a  cheby1(4..025.0.0394); 

postfilsig  a  filter(b2,a2,lim_sig);  %  Postfilterir)g  the  signal 

clear  llm_sig 

num  a  length(postfilsig); 

postfilslga[postfilsig(1,dd(j):num)  postfilsig(1,1:(dd(j)-1))];%  Filter 
sig_in  a  reshape(postfilsig,bit(j),m+1);  %  delay 

clear  postfilsig 
8ig_in  a  conj(sig_in’); 

rec_sig  a  demod(sig_in,m);  %  BPSK  demodulation 
dear  sig_in 

errorsO’) »  compare(md_o1(1:m-2),rec_sig(1:m-2));  %  Checking  errors 
clear  rec.sig 
end 

diary  juan.d 

errors  %  Saving  the  results  in  a  diary  file 
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diafy  off 

DBPSPO.M  (DBPSK  ADJACENT  CHANNEL  INTERFERENCE  POWER  COMPUTATION) 

%  This  M_rile  computes  the  power  in  the  main  channel  and  the  power  of  the 
%  adjacent  channel  that  is  leaking  into  the  main  channel.  After  that  a  ratio 
%  between  both  powers  is  obtained. 

m  s  1000;  %  Number  of  bits 

md_o  3  msg(l0,m);  %  Raixfom  message  generation 

md_o1  =  msg(25,m); 

dif_o  s  dif_cod(md.o);  %  Differentially  encoding  the  message 

dif_o1  =  dif_cod(md_o1); 

map_o  s  map(dif_o);  %  Mapping  the  message 

map_o1  =  map(dif_o1); 

delta_f  =  100000;  %  Separation  between  channels 
f1  =  0; 

clear  rrKl_o  dif_o 

t  =  1/384000;  %  Sampling  inten/al 

T  »  [(1/2400)  (1/4800)  (1/9600)  (1/19200)];  %  Bit  durations 

forj*1.4, 

mod_sig  =  modul(map_o,T(j),t.delta_f);  %  BPSK  modulation 

mod_sig1  =  modul(map_o1 ,  T(j),t.f1); 

mod_sig  *  mod_sig’; 

mod_sig  «  mod_sig(:); 

mod_sig  =  mod_sig’; 

mod_sig1  =mod_sig1'; 

mod_sig1  =  mod_sig1(:); 

mod_sig1  =mod_sig1’; 

[bl.al]  *  cheby1(6..01, 0.0651); 

prefiLsig  *  filter(b1  ,a1  ,mod_sig);  %  Prefiltering  the  signal 

prefiLsigl  =  filter(b1,a1,mod_sig1); 

lim.sig  s  limiter(prefiLsig);  %  Hard  limiting  the  signal 

Itfn_sig1  =  limiter(prefiLsigl); 

clear  prefiLsig  prefiLsigl 

[b2.a2]  =  cheby 1(4, .025, 0.0394); 

postfiLsig  =  fitter(b2,a2,lim_sig);  Postfiitering  the  signal 

postfiLsigI  =  fiRer(b2,a2,lim_sig1); 

clear  lim.sig  Iim_sig1 

II  s  length(postfiLsig); 

power_ad(j)  =  sum(abs(postfiL8>g)A2)/ll;  %  Computing  the  power 
power_base(j) «  sum(abs(postfiLsig1).^2)/ll; 
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dear  postfiLsig  postfiLsigI 
clear  mod.sig  mod.sigl 

r)orTn_power0  »  power_base(i)/|power_adO'):  %  Computing  the  ratio 
noiTn_pow_dB(j)  =  10*log10(norm_power(j)); 
end 

diary  juan.d 

norm_pow_dB  %  Saving  the  results  in  a  diary  file 
diary  off 


DEMOD.M  (DBPSK  DEMODULATION  FUNCTION) _ 

function  out  =  demod(in,m) 

%  This  M_file  performs  noncoherent  demodulation  of  DBPSK.  The  matrix  in 
%  contains  the  sampled  DBPSK  waveform  and  m  is  the  number  of  bits  that  this 
%  waveform  represents. 

o  =  ones(1,m): 
for  i3s2:m+1 , 

dif(i-l)  =  abs(sum(in(i,:))  -  sum(jn(i-1,:))); 
au(l-l)  =  abs(sum(in(i,;))  +  sum(ln(i-1 ,:))); 
metric(i-l)  *  dif(i-1)  -  su(i-l): 
if  metric(l-l)  <  0, 
o(i-l)  =  0: 
end 
end 

out  =  o; 


DIF_COD.M  (DIFFERENTIALLY  ENCODING  FUNCTION) _ 

function  dlf_o  =  dif_cod(ln) 

%  This  M_File  differentially  encodes  a  bit  stream  that  is  input  in  the 
%  variable  in. 

a  s  length(in); 
y  =  I1.2eros(1,a)l; 
for  isi.a, 

y(i+1)  *  xor(in(i),y(i)); 
erxf 

dif_o  =  y; 
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ESTBPS.M  (ESTIMATION  OF  CARRIER  TO  INTERFERENCE  RATIO] 

%  This  is  the  main  program  to  estimate  the  carrier  to  interference  ratio.  In 
%  this  case  the  modulation  used  is  DBPSK  but  the  method  holds  for  any 
%  modulation  scheme. 

m  s  1000;  %  Number  of  bits 

md_o  =  msg(40,m);  %  creating  the  random  message. 

itkI_o1  «  msg(6S,m);  %  creating  the  interfering  message. 

dif_o  s  dif_cod(md_o);  %differentially  erKX)ding  the  message 

dif_o1  ®  dif_cod(md_o1): 

map_o  s  map(dif_o);  %  mapping  the  message. 

map_o1  SB  map(dif_o1); 

T  =  [(1/2400)  (1/4800)  (1/9600)  (1/19200)1;  %  Bit  durations 
t  =  1/384000;  %  Sampling  interval. 
delta_f  s  100000;  %  frequency  separation 
f1=0; 

clear  md_o  dif_o  md_o1  dif_o1 
forjsBl;4, 

mod.sig  »  modul(map_o,T(j),t.delta_f);  %  OBPSK  modulation 

mod_sig1  =modul(map_o1.T(j),t,f1); 

inpass_b  bs  mod.sig  +  mod^sigl ;  %  Adding  both  messages 

inpass_b  =  inpass.b'; 

inpass_b  bs  inpass_b(;); 

inpass.b  =  inpass_b'; 

mod_slg1  simod_sig1'; 

mod_sig1  =  mod_sig1(:); 

mod_sig1  b:  mod_sig1'; 

[b1,a1]  *  cheby1(6,.01,0.0651); 

prefiLsig  « fiiter(b1  ,a1  .inpass.b);  %  Ritervig  both  messages 
prefii.sigi  *5  filter(b1,a1,mod_sig1);  %  Filtering  the  main  message 
lim.sig  bs  iimiter(prefiLsig);  %  Hard  limiting  both  messages 
clear  prefiLsig 

[b2,a2]  B>  chebyl  (4, .025, 0.0394); 
postfiLsig  sb  fiKer(b2,a2.lim_s{g); 
postfil_8ig1  sBfilter(b2,a2,preriLsig1); 
clear  lim.sig  prefiLsigI 
dd  s  sum(postfiL8ig1.^2); 

gainO) «  8um(postfil_slg.*postfll_8lg1)/dd;  %  computing  the  GAIN 


43 


in_estimate  »  gain(j)*postfiLsig1 ;  %  computing  baseband  estimate 
band_estimate  =  postfiLsig  •  in_estimate:  %  computing  the  interfering 
%  estimate 

clear  postfiLsig  postfiLsigI 
ff  s  ler)gth(in_estimate); 
power_in(j) » (sum(abs(in_estimate).^2))/ff 
power_band(j)  =  (sum(abs(barKl_estimate).^2))/ff 
C.toJO)  =10*iog10(power_in(j)/power_barKi(i));  %  computing  the  C/I 
clear  in_estimate  barxl.estimate 
clear  mod_sig  mod_sig1  inpass.b 
end 

diary  juan.d 

C_to J  %  Saving  the  results  in  a  diary  file, 
diary  off 


L1MITER.M  (HARD  LIMITER  FUNCTION) _ 

function  out=limiter(in) 

%  This  M_File  performs  a  hard  limiting  effect  over  a  modulated  signal.  This 
%  signal  is  contained  in  the  vector  in 
ss  s  abs(in); 
out  s  inVss; 


MAP.M  (MAPPING  FUNCTION) _ 

function  out  s  map(in) 

%  This  M_File  maps  a  bit  stream  to  0  or  pi  to  be  able  to  perform  afterwards 
%  a  BPSK  modulation. 

a  =  length(in); 
fori=1:a, 

If  (in(l) «  0), 
out(i)  *  0; 
else 

out(i)  =  pi; 
erxJ 
end 
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MODUL.M  (BINARY  PHASE  SHIFT  KEYING  MODULATION  FUNCTION] 

function  out »  modul(in,T,t,fc) 

%  This  M_File  pe  forms  BPSK  modulatioa.lt  accepts  the  signal  in,  the  bit 
%  duration  T,  the  sampling  interval  t  and  the  carrier  frequency  fc  as  inputs 
time  =  0:t:(T-t); 
a  s  length(in); 

P*  1: 
fors=1;a, 

time  =  time  +  (p  -  1)*T; 

P  =  2: 

if  in(s) »  pi, 

out(s,:)  =  exp(j*(2*pi*fc*time  +  pi)); 
else 

out(s,:)  =  exp(j*(2*pi*fc*time)); 
end 
er)d 

MSG.M  (MESSAGE  GENERATION  FUNCTION) _ 

function  u  -  msg(seed,k) 

%  This  M-file  accepts  a  data  vector  with  seed  for  rand  and 
%k  the  number  of  bits  that  will  be  returned  in  the  vector  u 
randfuniform’) 
rand(’seed',seed) 
u  =  round('-and(1,k)): 


XOR.M  (EXCLUSIVE-OR  FUNCTION) _ 

function  a  =  xor(in,in1) 

%  This  M_File  performs  the  xor  logic  operation 
if  in  Ks  ini, 
a  s  0; 

else 

a  ss  1; 

end 


45 


APPENDIX  B. 


CODEMOD.M  (COHERENT  MSK  DEMODULATION  FUNCTION) _ 

function  out  s  codemod(jn,t.T.fc,h,m) 

%  This  function  performs  coherent  demodulation  of  Minimum  Shift  Keying  using 
%  correlation,  sampling  and  integration  in  each  of  the  two  branches  of  the 
%  receiver.  The  integration  is  performed  over  a  period  equal  to  twice  the  bit 
%  duration  and  the  decision  is  made  by  alternatively  evaluate  the  output  of 
%  the  two  branches. 

time  »  0:t:(T-t); 
dd  =  1; 
ff  =  1. 
for  8=1  :m. 

phil  =  cos(pi*time*lVT); 
phi2  =  sin(pi*time*h/T); 
if  rem(s,2)  -=  0, 
vec1(cld,;)  =  in(s,:).*phi1; 
vec2(dd,:)  =  in(s,:).*phi2; 
dd  =  dd-«- 1; 
else 

vec3(ff,:)  =  in(s,:).*phi1; 
vec4(ff,:)  =  in(s,:).*phi2; 
ff  =  ff+  1; 
end 

time  =  time  +  T; 
erKl 

vec2  =  vec2  vec4; 
vec2  =  vec2’: 
sec  =  sum(vec2); 
last  =  8um(vec3(ff-1 
[rr  cc]  =  size(vecl); 
veci  =  vec1(2:rr,;): 
vec3  =  vec3(1:rr-1,;); 
veci  =  veci  +  vec3: 
vecl  =  veci’; 
one  =  sum(vecl); 
one  s  [one  last]; 
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H  real(one(ee))  >  0, 
estl(ee)  >0; 
else 

estl(ee)  spi; 
end 

M  jmag(sec(ee))  >  0, 
est2(ee) « -pi/2; 
else 

e$t2(ee) «  pi/2: 
end 
end 

if  est2(1)  ss -pi/2, 
dec(1)s0; 
else 

dec(1)  s  1; 
end 
k=1; 
v  =  1; 

for  ggs2:m, 
if  rem(gg-1,2)  -=  0, 

if  (est1(k)»0  &  est2(k)ss<pi/2)  I  (est1(k)»spi  &  est2(k)«spi/2), 
dec(gg)  =  0; 
else 

dec(gg)  =  1: 
erxl 

k  =  k+  1; 
end 

if  rem(gg-1,2)==0, 

if  (est1(v)«0  &  est2(v+1)==-pi/2)  I  (est1(v)=*pl  &  est2(v-«-1)==pi/2), 
doc(gg)  *  0; 
else 

dec(gg)*1: 

end 

v  =  v+  1; 
end 
end 

outs  dec; 
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CPFSKMOD.M  (CONTINUOUS  PHASE  FREQUENCY  SHIFT  KEYING  MODULATION 
FUNCTION)  _ 


function  out  =  cpfsKmod(in,T,t.fc.h) 

%  This  M_file  performs  the  modulation  of  CPFSK  with  any  modulation  index 
%  since  it  accepts  h  as  an  input. 
tetaO  s  0; 
a  s  length(in); 
time  =  0:t:(T -t); 
for  s=1  :a, 
if  s  ==  1 , 
teta  =  0; 
else 

tetaO  s  tetaO  +  in(s-l); 
teta  =  pi*h*teta0; 
end 

time  s  time  +  T; 
if  in(8)  ==  -1, 
f1  =fc-(h/(2*T)): 

mod_output(s,:)  =  exp(j*(2*pi*f1*time  +  teta  +  s*pi*h)); 
else 

f2  =  fc  +  (h/(2*T)); 

mod_output(s,:)  =  exp(i*(2*pi‘f2*time  +  teta  •  s*pi*h)); 
end 
end 

out  =  mod_output’; 
out  =  out(:); 
out  =  out’; 

EUCDIS.M  (EUCLIDEAN  DISTANCE  FUNCTION) _ 

function  D  =  eucdis(q,R) 

%  This  M-file  finds  Euclidean  distance  of  elements  in  vector  R  from 
%  q  unit  amplitude  vectors  equally  spaced  on  tue  unit  circle.  It  stores 
%  these  as  rows  of  D. 

L  =  length(R); 
nxiex  =  1  ;q; 
dph  =  2*pi/q; 

MO  s  exp(j*(dph.*(index-1))): 
for  1=1  :L, 
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D(l,:)  -  abs(R(0.'‘ones(MO)  -  MO); 


MAPPER.M  (MAPPING  FUNCTION) _ 

function  output «  mapper(in); 

%  This  M_File  maps  the  bit  stream  to  I's  or  >rs  to  be  able  to  perform  MSK 
%  modulation  afterwards, 
k  s  in  ss  0; 
k»  -k; 

output »  k  +  in; 


MATCH.M  (OFFSET  FUNCTION) _ 

function  [out.outi]  s  match(N,in,in1) 

%  This  M_File  matches  vectors  in  and  ini  which  are  offset  by  N  positions, 
if  length(in)  --  length(inl), 
out  a  in(1:iength(in)  •  N); 
out1  a  in1(N+1:length(in1)); 
end 


MSKVI.M  (MSK  RECEIVER  WITH  VITERBl  ALGORITHM) _ 

%  receiver  for  MSK  with  Viterbi  decoding 
clear 
m  a  1020; 
diary  juan.d 

md_o1  a  msg(40,m);  %  Creating  the  mesage 
map_o1  a  mapper(md_o1);  %  Maji^ir^g  the  function 
T  a  [(1/2400)  (1/4800)  (1/9600)  (1/19200)];  %  Bit  duration 
t  a  1/384000;  %  Sampling  intenral 
flaO; 

h  a  0.5;  %  Modulation  index 

sigma  a  [2*sqtt(2)  2  sqrt(2)  1];  %  StarKiard  deviation  of  the  noise 
clear  rrxf_o2  md_o3 
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forkks1:4, 

mod_sig1  =  cpfslonod(map_o1.T(kk),t,f1,h);  %  MSK  modulation 
ch_sig  s  awgn(mod_sig1  ,sigma(kk));  %  Adding  the  Gaussian  noise 

clear  mod_sig1 

dem.sig  >  videmod1(ch_sig,t.T(id(),h);  %  Mapping  the  signal  to  the 

%  euclidean  plarie 

clear  ch_sig 

TT  s  zeros(2,60);  %  Input  matrix  to  the  Viterbi  Algorithm 

vipath  s  [1  0  1  2  1  4;  1  1  2  2  0  3]: 

forqqs1:m, 

D  =  dem_alg(qq,;): 

TT  =  softv(1 ,2, 20, TT, vipath, D):  %  Viterbi  algorithm  function 
vLsig(qq)  =  TT(1,60): 
end 

clear  dem_sig 

[mes_o1  .rec_sigl  =  match(19.md_o1.vLsig):  %  Offset  function 
errors(kk)  =  compare(mes_o1,rec_sig);  %  Checking  enrors 
clear  rec.sig  mes_o1  vLsig 
end 

errors  %  Saving  the  results  in  a  diary  file 
diary  off 


REMSK.M  (COHERENT  MSK  MAIN  PROGRAM) _ 

%  receiver  for  MSK 
m  =  1000; 
diary  juan.d 

md_o1  =  msg(40,m);  %  Creating  the  message 
md_o2  =  msg(43,m);  %  Creating  the  interfering  message 
md_o3  =  msg(65,m); 

map_o1  s  mapper(md_o1);  %  Mapping  the  message 
map_o2  s  mapper(md_o2);  %  Mappkrg  the  interfering  message 
map_o3  s  mapper(md_o3); 

T  » [(1/2400)  (1/4800)  (1/9600)  (1/19200)];  %  Bit  durations 
t  s  1/384000;  %  Sampling  interval 
f1  =0; 

delta.f  =  100000;  %  Frequency  separation 
h  =  0.5;  %  Modulation  index 
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sigma  « [2*8qrt(2)  2  sqft(2)  1];  %  Standaid  deviation  of  the  noise 

bit -[160  80  40  20]; 

dd  -  [38  38  3841];  %  Filter  delay 

dear  md_o2  md_o3 

for  kk«1:4, 

mod_sig1 «  Cf>fskmod(map_o1,T(kk),t,f1,h);  %  MSK  message 
modulation 

mod_sig2  >  cpf8kmod(map_o2,T(kk),t,ddtaJ,h);  %  MSK  interfering 

message 

mod_sig3  >  cpf8i(rTK)d(map_o3.T(kk),t,-deltaLf.h);  %  modulation 
nfKxf.sig  >  mod.sigl  +  mod_8ig2  +  mod_8ig3;  %  Adding  the 
messages 

dear  mod.sigi  mod_8ig2  mod_sig3 

ch_8ig  «  awgn(mod_sig,sigma(kk));  %  Adding  the  Gaussian  noise 

dear  rrKxf.sig 
ch_8ig  =  ch_aig’; 
ch_slg  «  ch_8ig(:); 
ch_sig  «  ch_sig'; 

[b1,a1]  =  cheby1(8..01 ,0.0851); 

prefil_sig  « filter(b1  ,a1  ,ch_sig);  %  PrefiRervig  the  signal 

dearch_sig 

Hm.sig  « limiter(prefiLsig);  %  Hard  limiter  effed 
dear  prefiLsig 

[b2.a2]  »  chebyl  (4, .025, 0.0394); 

postfiisig  B  fiiter(b2,a2,lim_sig);  %  Postfiltering  the  signal 

dear  lim_sig 

num  « length(postfilsig); 

postfiisig  B  [postfilsig(1,dd(j):num)  postr)isig(1,1:dd(j)*1)];  %  Filter 
sig_in  s  reshape(postfilsig,bit(kk),m);  %  delay 

dear  postfiisig 
sig_in  »  conj(sig_in’); 

rec_sig  b  codemod(sig_in,t,T(kk),f1,h,m);  %  Coherent  demodulation  of 

%  MSK 

dearsig_in 

errors(kk)  b  oompare(nxi_o1.rec_8ig);  %  Cheddng  errors 
dear  rac_sig 
erKf 

errors  %  Savirni  the  result  in  a  diary  file 
diary  off 
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SOFTV.M  fSOFT  VITERBI  DECODER] 


function  PHN  *  softv(k,K,Np.PH.T.D) 

%  Soft  Viterbi  Decoder 

%  Paul  H.  Moose 

Univ.  degli  Studi  di  Padova 
%  17-05-91 

%  This  M-file  decodes  k  bit  msgwords  from  2^  real  metrics 
%  (These  may,  for  example,  represent  the  "distance"  of  the 
%  received  modulation  value  from  each  of  2^  modulation 
%  values.) 

%  The  state  transition  information  for  a  2'^K  state  trellis  is  in 
%  the  2^K  by  3*2'^k  matrix  T.  Each  of  the  2'^k  entering  paths  to 
%  each  state  has  its  source  state  (one  of  2^K),  path  msgwords  (one 
%  of  2'^k)  and  path  codeword  (one  of  2^)  listed  in  the  state  row. 

%  The  path  histories  are  kept  in  matrix  PH  that  is  2''K  by  3*Np. 

%  The  path  history  for  each  state  contains  source  state,  path 
%  weight  and  path  codeword  for  Np  previous  states. 

%  The  output  PHN  is  the  update  of  PH,  the  new  path  history. 

%  The  decoded  codeword  is  in  the  last  column  of  PHN.  (They  should 
%  "merge"). 

%  The  past  histories  are  updated  on  the  basis  of  the  "minimum 
%  metric".  You  can  change  this  to  the  "maximum  metric"  if  desired  as 
%  indicated  in  the  comments  in  the  code. 

% 

forffsV.K 

X(ff.2)  *  D(T(ff,3))  +  PH(T(ff,1).2);  %path  weight 
X(ff,1)  =  %path  source  state 

X(ff.3)  =  T(ff,2):  %path  code  word  T(ff,3).Chg  to  T(ff.2)  for  msgword 
for  1=2:2^ 

wt  =  D(T(ff,3*l))  +  PH(T(ff.3*l-2).2); 
if  wt  <  X(ff,2)  %  The  <  selects  mn  metric 
X(ff.2)  =  wt; 

X(ff,1)  =  T(ff, 3*1-2); 

X(ff,3)  *  T(ff.3*H);  %  Ghg  to  T(ff,3*l)  for  codeword 
end 
end 

%  We  need  now  to  append  old  paths  to  new  paths  to  get  survivors. 

PHN(ff.:) » [X(ff.:)  PH(X(ff.1),1:3*Np-3)l; 
er>d 


VIDEMOD1.M  (VITERBI  DEMODULATION  FUNCTION) 

function  out »  videmod1(in,t,T,h) 

%  This  M_File  accepts  a  modulated  signal  and  matches  it  on  the  euclidean 
%  plane.  The  euclidean  distance  from  these  points  to  4  different  points 
%  is  found  and  the  metric  is  returned  to  be  used  as  an  input  in  the  soft 
%  Viterbi  decoder 

no_mat «  Q; 
map  » [-1  1]: 

[n  cc]  s  size(in); 
time  »  0:t;T-t: 
for  sssi  :rr, 

in(ss,:) « in{ss,:)*exp(j*(ss-1)*pi*h); 
for  m_ary=1 :2, 

XX  ss  exp(j*pi*h*time*map(m_ary)/T); 
first  =  xx*conj(in(ss,:)'): 
no_mat  *  Ino_mat  first]; 
end 

demod(ss,:)  =  real(rK>_mat); 
no_mat  =  Q; 
erxi 

R  *  j*d0mod(:,1)  +  demod(;,2); 
out «  eucdis(4,R); 
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APPENDIX  C. 


CPFSDIS.M  (EUCLIDEAN  DISTANCE  FUNCTION  FOR  CPFSK  WITH  h  =  0.4) 
function  D  «  cpfsd»(R) 

%  This  M-file  finds  Euclidean  distance  of  elements  in  vector  R  from 
%  10  unit  amplitude  vectors  on  the  unit  circle  (This  is  the  case  of  CPFSK 
%  with  h=0.4).  It  stores  these  as  rows  of  D. 

L  s  length(R); 

dph  s  [0.235  1.336  1.183  2.672  2.867  -2.221  -2.491  -1.296  -1.101 
0.388]; 

MO  s  exp(j*dph); 
forl»1:L, 

0(1.;)  >  abs(R(l).‘ones(MO)  -  MO); 
end 


CPFSK.M  (VITERBI  ALGORITHM  RECEPTION  OF  CPFSK) _ 

%  receiver  for  CPFSK  with  h=  0.4 
m  = 1000; 
diary  juan.d 

md_o1  =  msg(40,m);  %  Creatbig  the  message 
map_o1  =  mapper(md_o1);  %  Mapping  the  message 
T  =  [(1/2400)  (1/4800)  (1/9600)  (1/1^00)];  %  Bit  durations 
t  s  1/384000;  %  Sampling  interval 
f1  =0; 

h  =  0.4;  %  Modulation  index 

sigma  s  [2*sqrt(2)  2  sqrt(2)  1];  %  Standard  deviation  of  the  noise 

clear  md_o2  md_o3 

forkk=1:4, 

mod_sig1  =  cpfskmod(map_o1,T(kk),t,f1,h);  %  CPFSK  modulation 
ch_sig  =  aw^(mod_sig1,sigma(kk));  %  Adding  the  Gaussian  noise 

clear  mod_sig1 

dem.sig  s  videmod(ch_sig,t,T(kk),h);  %  Mapping  the  signal  in  the 

%  euclidean  plane 
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dear  ch_8ig 

TT  M  zeros(5,60);  %  Input  matrix  to  the  Viterti  algorithm 
vipath-{2  0 3 5  1  10;  3  05  1  1  2;  4  072  1  4;  5  09  3  1  6;  1  0  1  4  1  8]; 
for  qq*1  :m, 

D  ■  dom_8ig(qq.:): 

TT  a  8oftv(1 ,5,20,TT,vipath,D):  %  VitMbi  decodng  function 
vi_8ig(qq)aTT(1.60); 
end 

dear  dem_8ig 

[me8_o1,rac_aig]  =  match(19,md_o1,vi_8ig):  %  Offaet  function 
erTors(kk)  =  compare(mea_o1  ,rec_8ig);  %  Checking  errors 
clear  rec.sig  mes_o1  vi_sig 
erKi 

errors  %  Saving  the  results  in  a  diary  file 
diary  off 


VIDEMOD.M  (VITERBI  DEMODULATION  FUNCTION) _ 

fundion  out  =  videmod(in,tJ,h) 

%  This  M_File  accepts  a  modulated  signal  and  matches  it  on  the  euclidean 
%  plane.  The  euclidean  distance  from  these  points  to  10  different  points 
%  is  found  and  the  metric  is  returned  to  be  used  as  an  input  in  the  soft 
%  Viterbi  decoder 

rx)_mat  *  Q; 
map  a  [-1  1}; 

[rr  cc]  a  size(in); 
time  a  0:t:T-t; 
for  ssal  :rr, 
for  m_aryai  .2, 

XX  a  oxp(j*pi*h*time*map(m_ary)/r): 
first  a  xx*coni(in(ss,:)’); 
rK)_mat  a  [no_mat  first]; 
erxl 

demodfss,;)  a  real(r)o_mat); 
no_mat  a  Q; 
end 

R  a  j*demod(:,1)  +  demod(:,2); 
out  a  cpfsdis(R); 
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APPENDIX 


D. 


AND.M  (AND  GATE  FUNCTION) _ 

function  out  =  and(in,in1) 

%  This  M  file  accepts  two  bits  as  inputs  and  perfonns  the  logical  "and"  operation 
%  between  them. 

if  in  ss  1  &  ini  1 , 

out  s  1; 
else 
out  =  0; 
end 


PECISION.M  (DECISION  BLOCK  FOR  NONCOHERENT  MINIMUM  SHIFT  KEYING) 

function  outsdecision(in) 

%  This  M_file  accepts  a  vector  that  represents  the  output  of  the  integrator 
%  in  noncoherent  reception  of  MSK  and  decides  whether  this  output  corresponds 
%  to  a  zero  or  a  one, 

a  =  length(in); 
out  =  zeros(1  ,a); 
for  j=1  :a, 
if  in(j)  >  0, 
out(j)  =  1; 
end 
end 


MSKDEMOD.M  (NONCOHERENT  MSK  DEMODULATION  FUNCTION) _ 

function  out  =  mskdemod(h) 

%  This  M_File  performs  noncoherent  MSK  demodulation  over  a  signal  contained 
%  in  the  matrix  in 

for  h=2:length(in). 
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y(h.;)  -  real(in(h.:).*coc^in(h-1.:)).*exp(-j*(pi/2))): 
end 
out »  y; 


PARITY.M  (PARITY  BIT  FUNCTIONt _ 

function  outi  »  parity(in) 

%  This  M_File  obtains  a  parity  bit  from  a  MSK  signal.  This  parity  bit  is 
%  going  to  be  used  in  the  single  error  correction  circuit, 
ss  s  length(in); 
for  hs3;ss, 

y1(h.;) » real(in(h.:).*  conj(in((h-2),;))): 
end 

outi  =  y1 : 


RENCMSK.M  (NONCOHERENT  RECEPTION  OF  MSK) _ 

%  receiver  for  NON  COHERENT  MSK 
diary  juan.d 
m  s  1005; 

md_o1  =  msg(40,m);  %  Creating  the  mah  message 
nid_o2  s  msg(43,m);  %  Creating  the  interfering  messages 
md_o3  =  msg(65,m); 

map_o1  s  mapper(md_o1);  %  Mapping  the  message 

map_o2  s  mapper(md_o2);  %  Mapping  the  interfering  messages 

map_o3  =  mapper(rTKl_o3): 

dd  s  [71  55  45  43];  %  Filter  delays 

T  «  [(1/2400)  (1/4800)  (1/9600)  (1/19200)];  %  Bit  durations 

t  =  1/384000;  %  Sampling  inten/al 

f1  =0; 

deltaj  s  100000;  %  Frequency  separation 
bits  [160  80  40  20]; 

si^a  s  r2*sqtt(2)  2  sqrt(2)  1];  %  Starxlaid  deviation  of  the  noise 
h  s  0.5,  %  Modulation  irKlex 

B  s  [0.0138  0.0275  0.055  0.1 1];  %  Bandwidth  of  the  first  Butterworth 
filter 

B1  s  [0.0188  0.0375  0.075  0.15];  %  Bandwidth  of  the  sec.  Butterworth 
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%  filter 

clear  md_o2  md_o3 
for  j=1 :4, 

mod_sig1  «  cpfsknK)d(map_o1,T(j),Ul.h):  %  MSK  modulation 

mod_sig2  »  cpfsKniod(map_o2.T0)>td^-f.h); 

mod_sig3  ®  cpfskmod(map_o3,T(j).t.*<*®lteJ.h): 

mod.sig  s  mod.sigi  +  mod_sig2  +  mod.sigS;  %  Adding  the  signals 

clear  mod_sig1  mod_sig2  mod_sig3 

ch.sig  s  awgn(mod_sig,sigma(j)):  %  Adding  the  Ge  jssian  noise 

clear  mod.sig 

ch_sig  s  ch_sig’; 

ch_sig  =  ch_sig(:): 

ch_sig  =  ch_sig': 

P31  ,a1]  a  chebyl  (6, .01 ,0.0651 ); 

prefiLsig  =  filter(b1  ,a1  ,ch_sig);  %  Prefilteimg  the  signal 

clear  ch_sig 

lim_sig  =  iHniter(prefiLsig);  %  Hard  limiter  effect 
clear  prefiLsig 

[b2,a2]  s  chebyl  (4,0.025,0.0394); 

postfilsig  =  fiiter(b2,a2,lim_sig);  %  Postfiltering  the  sigruU 

clear  lim_sig 

[b3,a3]  s!  butter(4,B(j)); 

recfil  =  filter(b3,a3,postfilsig);  %  Butterworth  filter  in  the  receiver 

clear  b3a3  postfilsig 

shapesig  =  reshape(recfjl,bit(j,  m); 

shapesig  s  conj(shapesig'); 

demsig  s  mskdemod(shapesig);  %  MSK  demodulation 

parsig  s  paiity(shapesig);  %  Parity  bit  creation 

clear  recfil  shapesig 

demsig  =  demsig'; 

demsig  »  demsig(:); 

demsig  s  demsig'; 

parsig  s  parsig'; 

parsig  »  parsig(:); 

parsig  s  parsig'; 

[b4,a4]  -  butter(4,B1(j)); 

demsig  » fHter(b4,a4,demsig);  %  Second  Butterworth  filter 
parsig  s  filter(b4,a4,parsig); 

demsig  =  [demsig(1,dd0'):length(demsig))  demsig(1,1:dd(j)-1)];  %  Filter 
parsig  [parsig(1,dd(j):length(parsig))  parsig(1,1:dd0-1)];  %  ctolay 
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demsig  « re8hape<dem8ig,bit(j),m); 
parsig  » reshape(parsig.bitG).ni): 
shapedem  >  sitfn(dem8ig); 
shapedeml  »  8hapedem(3:m); 
shiq;>apar «  8Ufn(par8ig); 

8hapepar1  >  8hapepar(3:m); 
dear  derraig 

dataout  s  deci8ion(8hapedeni1);  %  Decision  block 
paraout «  decision(shapeparl); 

Sd  Single  error  correction  circuit 

datar  s  sinerror(dataout(2:length(dataout)),paraout(2:length(paraout))); 
clear  dataout  paraout 

clear  shapedem  shapedeml  shapepar  shapeparl 

errorsG) »  compare(md_o1(4:m-2),datar(2:length(datar)));  %Checking 

%  errors 

clear  datar 
end 
errors 

%  Saving  the  results  in  a  diary  file 
diary  off 


SINERROR.M  (SINGLE  ERROR  CORRECTION  FUNCTION) _ 

function  correct »  sinerror(in,in1) 

%  This  M_File  performs  a  single  error  correction  accepting  as  inputs  a  vector 
%  in  which  contains  the  data  and  a  vector  ini  which  contains  the  parity  bits 
in  -  fO  0  0  in]; 
ini  » [0  0  0  ini]; 
out  =  zeros(1  ,length(in)); 
out3  >  zeros(1  ,length(in)); 
out2  s  zerDS(1,length(in)); 
for  kks3:iength(in), 

out(kk)  s  xor(in(kk),in(kk-1)); 
outl(kk)  s  xor(out(kk),in1(kk)); 
out2(kk)  s  xor(out1(kk),out3(kk-1)); 
out3(kk) «  and(out2(kk-1),outi(kk)); 
correct(kk-1)  >  xor(out3(kk),in(Kk-1)); 
end 

correct »  correct(4:length(in)-1); 


59 


LIST  OF  REFERENCES 


1 .  Couch,  Leon  HI,  Digital  and  Analog  Communications  Systems,  McMillan,  1 990. 

2.  Ha,  Tri  T.,  Digital  Satellite  Communications,  McGraw  Hill,  1990. 

3.  Hughes  Aircraft  Company  Space  and  Communications  Group,  "Proposal  for  UHF 
Follow-On  Communications  Satellite,  Vol.  D,"  January,  1988. 

4.  Murota,  K.,  and  Hirade,  K.,  "GMSK  Modulation  for  Digital  Mobile  Radio 
Telephony,"  IEEE  Trans.  Commun.,  Com-29,  No.  7,  July,  1981. 

5.  Kahnaman,  Betulham,  "Performance  Evaluation  of  UHF  Fading  Satellite  Channel 
by  Simulation  for  Different  Modulation  Schemes,"  Master's  Thesis,  Naval 
Postgraduate  School,  Monterey,  CA,  December,  1992. 

6.  Haykin,  S.,  Communications  Systems,  2nd  ed.,  Wiley  and  Sons,  1983. 

7.  Proakis,  J.G.,  Digital  Communications,  2nd  ed.,  McGraw  Hill,  1989. 

8.  Mason,  Lloyd  J.,  "Analysis  and  Simulation  of  a  DMSK  Receiver,"  SCC 1983 
Conf.  Rec.,  Ottawa,  June,  1982. 

9.  Masamura,  T.,  Samejima,  S.,  Morihiro,  Y.,  and  Fuketa,  H.,  "Differential  Detection 
of  MSK  with  Nonredundant  Error  Correction,"  IEEE  Transactions  on 
Communications,  Vol.  Com-27,  No.  6,  June,  1979. 

10.  Bhargava,  Haccoun,  Matyas,  and  Nuspl,  Digital  Communications  by  Satellite,  1st 
ed.,  John  Wiley  and  Sons,  1981. 

11.  Crozier,  S.,  Mazur,  B.,  and  Matyas,  R.,  "Performance  Evaluation  of  Differential 
Detection  of  MSK,"  GLOBECOM 82  IEEE  Telecommun.  Conf.  Rec.,  Miami, 
Nov.-Dee.,  1982. 


12.  deBuda,  Rudi,  "Coherent  Demodulation  of  Frequency-Shift-Keying  with  Low 
Deviation  Ratio,"  IEEE  Transactions  on  Communications,  June,  1972. 


60 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Infonnation  Center 
Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library,  Code  52 

Naval  Postgraduate  School 
Monterey,  CA  9394J-5000 

3 .  Director  Space  and  Electronic  Combat  Division  (N64) 
Space  and  Electronic  Warfare  Directorate 

Chief  of  Naval  Operations 
Washington,  DC  93943-5000 

4.  Chairman,  Electronic  Warfare  Academic  Group,  Code  EW 
Naval  Postgraduate  School 

Monterey,  CA  93943-9528 

5 .  Professor  P  H.  Moose,  Code  EC/Me 
Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

6.  Professor  R.  Clark  Robertson,  Code  EC/Rc 
Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

7.  Servicio  De  Guerra  Electronica 
Comando  De  Operaciones  Navales 
Puerto  Belgrano,  CP  8000,  Argentina 

8.  Sevicio  De  Inteligencia  Naval 
Estado  Mayor  General  De  La  Armada 
Buenos  Aires,  Argentina 


61 


9.  Direccion  De  Instruccion  Naval 
Estado  Mayor  General  De  La  Armada 
Buenos  Aires,  Argentina 


t 


